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dispersion  curves  for  this  plate.  The  flexural  (acoustic)  branch  has 
phase  and  group  velocities  approaching  the  shear  velocity  of  the  plate 
material  as  frequency  Increases.  Stuart  [12-14]  explored  the  solution 
for  the  same  plate  further  with  new  Insights  into  the  leaky  waves 
emanating  from  the  plate  and  obtained  a  more  accurate  solution  when  the 
angle  of  observation  approaches  the  coincidence  angle.  This  new  solution 
would  be  accurate  at  closer  observer  distances  than  Felt's  [11]. 

Sound  radiation  from  beam-reinforced  plates  excited  by  point  or 
line  forces  has  been  investigated  by  few  authors.  Romanov  [15]  obtained 
the  solution  of  the  radiated  pressure  from  a  plate  reinforced  with  beams 
and  excited  by  a  line  force.  Felt  and  Saurenman  lib]  analyzed  the 
acoustic  radiation  of  a  point-excited  plate  reinforced  with  beam,  but 
confined  their  interest  to  high  frequencies.  Gorman  [17]  obtained  the 
solution  for  a  plate  reinforced  by  many  beams  and  excited  by  a  line  force 
parallel  to  beams.  Ills  solution  thus  makes  the  beam's  influence  purely 
as  rotatory  and  transverse  impedances  with  no  waves  travelling  in  the 
beams.  Garrellck  and  Lin  [18]  analyzed  the  radiation  from  beam- 
reinforced  plate  and  confined  their  attention  to  on-axis  response.  Lin 
and  Hayek  [19]  obtained  the  exact  solution  for  the  radiation  from  a 
po int- exc i ted  plate  reinforced  by  one  beam. 

Investigations  into  the  transient  acoustic  radiation  from  a 
plate  under  the  Influence  of  time  dependent  point  forces  were  investi¬ 
gated  by  Magrab  and  Reader  [20]  and  Stuart  [12].  Magrab  and  Reader  |20] 
predicted  the  tin*'  signature  of  the  radiated  pressure  for  a  CW  pulse. 
However,  the  solution  was  valid  only  after  the  acoustic  arrival  and 
there  was  an  error  in  the  choice  of  the  poles  in  the  solution.  Stuart 
[12]  predicted  the  Impulse  response  of  an  elastic  plate.  His  formulation 


accounts  tor  shear  anil  rotarv  Inert  l.i  effects  ot  the  plate,  .uni  the 
fluid  lo.ul  i  ue.  ooupline,  ot  the  acoustic  medium.  The  resu  1 1  lin;  far-tield 
radiated  pressure  .is  determined  from  the  standard  saddle  point  method  Is 
tor  t  imos  hotore  and  utter  the  aeons  tie  arrival.  However,  the  origin  ot 
the  tlrst  arrival  is  not  prodiet able  by  his  method.  The  solution  alter 
the  aeoust ie  arrival  was  predieted  to  he  mono t out ea 1 1 v  decaying  solut ion 
However,  it  is.  more  phvsieallv  reasonable  to  assume  that  the  pl.ite  will 
vibrate  treelv,  genet  at  ine,  an  ose  i  1  la  t  i  up,  and  decaying  time  signature. 

1  .  .'  (lea  I  s  a  1  t  he  S  t  ml  v 

lleeause  ot  the  inaeeurae  i  es  or  ambiguities  in  the  prodiet  ion  of 
the  time  signature  i>t  the  radiated  pressure  emitted  t rom  a  tiansientlv 
oxeited  pl.ite,  this  study  undertakes  to  prediet  the  solution  ot  stieh  a 
problem  exaetly.  Kurt hermore ,  sinee  the  exaet  solution  is  bound  t o  be 
romp  I  ieat ed .  approximate  solutions  will  be  sought  that  are  more  .leeurate 
than  the  em  rent  Iv  available  ones. 

1  ho  t  I  mi'  s  i  gnu  t  u  r  e  ot  tin*  radiated  pressure  is  t  o  bo  oomputed  tor 
an  impulse  toree,  a  constant  force,  a  Mijuare  pulse,  and  a  ilW  pulse.  The 
behavioi  ot  the  solutions  at  short  times  and  tor  lone,  times  will  be 
exam i nod . 


( 
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CHAPTER  II 


STATEMENT  OF  THE  PROBLEM  AND 
THE  MATHEMATICAL  MODEL 

2.1  St  atement  of  the*  Problem 

Consider  an  infinite  thin,  undamped  elastic  plate  excited  by  a 
point  force  as  shown  in  Figure  2.1.  The  plate  is  in  contact:  with  an 
infinite  acoustic  medium  occupying  the  space  z  >  0  and  has  vacuum  in 
the  space  z  <  0  .  The  point  force  is  applied  dynamically  at  t  =  0 
with  a  prescribed  amplitude  versus  time  history,  generating  transient 
elastic  waves  in  the  plate  and  transient  acoustic  pressure  in  the  fluid. 
The  transient  radiated  pressure  time  history  is  to  be  evaluated  at 
different  points  in  the  acoustic  medium. 

2 . 2  Matjiema  t  icn  1  Model 

An  infinite  thin  elastic  plate,  initially  undeformed  and  at  rest 
is  excited  by  a  dynamic  point  force  as  shown  in  Figures  2.1  and  2.2. 
Because  of  the  axisymmetry  of  the  force  field,  i.e.,  a  point  force,  the 
geometric  configuration  can  be  expressed  in  cylindrical  coordinates 
independent  of  circumferential  angle  0  .  The  assumption  of  thin 
lossless  elastic  plate  allows  one  to  use  the  Euler  plate  theory  as  the 
mathematical  model  for  the  governing  equations  of  motion.  However,  the 
assumption  of  a  thin. plate  breaks  down  when  the  wavelength  in  the  plate 
is  less  than  eight  times  the  thickness.  Thus,  the  Euler  plate  theory 
would  lead  to  incorrect  group  and  phase  velocities  for  wavelengths 
shorter  than  eight  times  the  thickness.  Therefore,  one  would  expect  a 
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slight  error  ii\  the  prediction  of  the  first  arrival  which  corresponds  to 
an  infinite  frequency  or  zero  wavelength.  The  equation  of  motion  for 
the  classical  plate  under  the  influence  of  a  force  field  can  he  written 


v'*U(r,t)  +  \  ~  (r.t)  =  -1-  If(r.t)  -  P(r,o,t)l  ,  (2.1) 

c  he  t 
P 


wlie  re 


•)  y 

.13 

-2  r  3r 
on 


t  >  0 
r  >  0 


where  W(r,t)  is  the  transverse  displacement  of  the  plate,  f(r,t)  is 
the  applied  force  field,  P(r,o,t)  is  surface  acoustic  pressure, 

P  =  (Eh  )/ll2(l  -  V“) ]  is  the  plate  flexural  rigidity,  c^  =  E/ [ 12  ( 1  — \r )  p  ] 
is  the  plate  velocity,  E  is  the  Young's  modulus,  h  is  the  thickness, 
p  is  the  plate  mass  density,  and  V  is  the  Poisson's  ratio. 

The  velocity  potential  <I>(r,z,t)  of  the  acoustic  medium  is 
governed  by  the  wave  equation: 


2  1  3“'<h 

V  >K r , z , t )  -  — j  — j  (r ,z,t)  =  0 

c  *■  3t 
o 


(2.2) 


7  7 

d*  13,  3 
3r2  r  K  3*2 


and  c  is  the  phase  velocity  of  the  acoustic  medium.  The  acoustic 
o 

pressure  P(r,z,t)  in  the  fluid  medium  is  related  to  the  velocity 
potential  by: 


P (r , z  ,  t )  = 


(2.3) 


1 
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where  0  is  the  mass  density  of  the  acoustic  medium.  The  boundary 

O 

condition  to  be  satisfied  is  continuity  of  normal  velocity  at  fluid- 
plate  interface,  i.e. , 


3/ 


( r  ,  o ,  t ) 


dW 

i)t 


(r.t) 


(2. 4) 


The  acoustic  medium  and  plate  are  considered  initially  undisturbed  and 
at  rest. 

In  addition  to  Equation  (2. A),  the  outward  radiation  condition 
must  be  satisfied. 


The  radiated  pressure  from  an  infinite  elastic  plate  excited  by 
an  impulsive  point  force  will  be  obtained  by  using  transforms.  The 
Fourier  complex  transform  is  applied  to  the  temporal  variable  t  ,  and 
the  Hankel  transform  to  the  spatial  variable  r  .  These  transform  pairs 
are  given  by  the  following  relations;  the  Fourier  complex  transform  pair 
is : 


* 

F  (u<) 


r  /  V  i<‘)t  . 

f  ( t )  e  d  t 


and 


f(t) 


1 

2tt 


F*((o)e"Uotd<J 


(2.5a) 

(2.5b) 


and  the  Hankel  transform  pair  is: 


F(k) 

and 

f(r) 


CV> 

r  f  ( r )  J  ( k  r )  d  r 
0 

o 


CO 

f  _ 

j  kF(k)  .l(>(kr)dk 
o 


(2. ha) 


(2 . bb) 
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2.3  Trans  formed  Equations 

Applying  Fourier  and  Hankel  transforms  to  Equations  (2.1),  (2.2), 
(2.3),  and  (2.4),  the  resulting  transformed  equations  are,  respectively: 


kSv*(k,w)  +  [-w2W  (k  ,oo)  ]  =  ^  [f  (k  ,oj)  -  P  (k,o,w)]  , 

h  c 

P  (2.7) 


2—* 

±4-  (k,z,u)  -  (k2  -  k2)  ?*(k,z,u)  =  0 

,  2  o 

dz 


(2.8) 


P* (k , z ,w) 


-iupo  <*>  (k,z,w) 


(2.9) 


(k,o,<a)  =  iwW 


(2.10) 


Equating  z  to  zero  in  Equation  (2.9),  and  substituting  into 
Equation  (2.7),  results  in  the  following  coupled  algebraic  equation: 


kV'(k.uj)  +  [-aj2w''(k,cj)]  =  ^  [f*(k,u)  +  iup  <J>*(k,o,U))] 

h  c 

p  (2.11) 


Solving  Equation  (2.11)  for  the  transformed  plate  displacement  W  (k,w), 
one  obtains: 


W  (k,w) 


- A  — & 

f  (k,w)  +  iU)p  $  (k,o,w) 
o 


(2.12) 


The  applied  force  field  to  be  considered  in  this  study  is  assumed  to  be 
a  point  force  acting  at  the  origin  of  the  coordinate  system.  Hence,  the 
applied  force  can  be  represented  by: 
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f(r.O  '  2(-^> 


(2.13) 


where  iS(r)  is  the  Dirac  delta  function,  and  Q(t)  represents  the  time 
dependent  amplitude  of  the  force.  Again  applying  Fourier  and  Hankel 
transforms  to  the  force  field  in  Equation  (2.13),  one  obtains: 


7*  ,  q(M>) 

t  (k  ,uj)  =  -V- 


(2.14) 


The  general  solution  to  Equation  (2.8)  can  be  written  as: 


-*  Vz  I’ 2 

i>  (k,z,.o)  =  Ae  +  Be 


>  0 


(2.15) 


where 


am! 


(k  -  k')x/~ 
o 


k  —  u)/  c 

o  o 

n  "St 

The  two  solutions  of  *P  (k,z,w)  represent  progressive  waves  that  travel 
outward  and  inward  in  the  z-direction.  Imposing  outward  radiation 
condition  on  Equation  (2.15)  yields: 


and 


-*  Fz 

(k  ,  z  ,u>)  =  Ae1 


2  1/2 

l  mi’  =  Im|k“  -  k  ]  >0 

o 


z  >  0 


(2.10) 


(2.17) 


Applying  Fourier  and  Hankel  transform  to  Equation  (2.4),  and  then 
equating  with  Equation  (2.16),  one  obtains: 


—*  AF 

W  (k,u>)  = 


(2.18) 
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Solving  for  A  from  Equations  (2.12)  and  (2.18),  one  obtains 
an  expression  for  the  pressure  amplitude  transform: 


V  2 

2TT 

1  4  cj 

_  “  P0 

"*  on 

+  ° 

hV 

l  PJ 

(2.19) 


From  Equations  (2.9),  (2.14),  (2.16),  and  (2.18),  the  transform  of  the 
pressure  field  is  given  by: 


P*(k,z,w)  = 


2  Tz 

P  u)  Q(ca)e 


!  4  w2  w2p  1 

2TTD,  (k“  -  -~)  r  + 

I  h  c 

P 


,  z  >  0  .  (2.20) 


CHAP  I  KK  I  I  1 


1 N  V  KK S  V.  X  PANS  FORMAT  1  ON 

3.1  Inverse  Fourier  Transform 

One  has  a  choice  of  which  of  the  inverse  transforms  to  perform 
first.  It  was  decided  to  perform  the  inverse  Fourier  transform  exactly 
first  so  that  the  integrand  in  the  llankel  inverse  transform  is  exact. 

Any  approx imations  performed  in  the  llankel  inverse  transformation  would 
thus  be  of  known  order.  However,  taking  the  llankel  inverse  transform 
first  and  performing  an  approximation  in  its  evaluations,  as  was  done  by 
Stuart  [12],  would  lead  to  an  approximate  integrand  in  the  Fourier 
inverse  t runs torn  and  hence  would  lead  to  an  unknown  order  of  approxima¬ 
tion.  Furthermore,  the  first  arrival  time  would  be  unknown  and  the 
solution  before  the  acoustic  arrival  would  not  have  an  initial  value. 

In  addition,  the  solution  after  the  acoustic  arrival  would  show  an 
exponentially  decaying  signature,  which  is  contrary  to  known  experi¬ 
mental  data. 

Applying,  inverse  Fourier  transform  on  time  to  liquation  (3.20), 
one  obtains: 

CO 

P  (k , ,  t )  =  2V  PV(k,.’.,w)e~UOtdu)  .  (3.1) 

4 

—<X) 


— CO 
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Using  the  convolution  theorem.  Equation  (3.2)  can  be  written  as: 


P(k , 2 , t ) 


where 


I*(k,z,t) 


t 

f  * 

j  I  (k,2,i )Q(t-T>di 
o 


co 


Tz-i^jt  , 
e  dcz 


2 


-0— )  r  + 

P 


^2Pp 

D 


(3.3) 


(3. A) 


is  the  impulse  response  of  the  system.  It  will  be  convenient  to 
normalize  the  angular  frequency  to  that  of  the  classical  coincidence 
frequency  defined  by: 


u 

o 


c 

c 


p 


such  that  the  normalized  frequency  £2  is  defined  as: 


P  = 


Thus,  Equation  (3.4)  can  be  rewritten  as: 

~  fz-itu  Pt 


I  (k, z ,  t ) 


8a  C>  0 


2 

47T“h 


P  e 


dP 


—CO 


(k4  -  a4p2)  r  +  a4P2$ 


(3.5) 


or,  by  dividing  the  integration  path.  Equation  (3.5)  becomes: 


I*(k,z,t) 


o  _ 

~  tz-iu)  Pt 
R"  e  °  dP 

(F4  -  a4p2)  7  +  a4p2B 

J 

— 


0  Tz-iui  Pt 
P“e  °  dP 


—4  4  2  4  > 

(k  -  a  P  )  T  +  a  P  S 


(3.5a) 


o 
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where 


k  =  kh  and  T  =  Fh 


Equation  (3.5a)  may  be  written  as: 
5 


I 


where 


2 

4rr“h 


) 


(3.6) 


and 


„  rz-io)  fit 
fi2  e  °  dfi 


—4  4  9  —  49 

(k  -  u  fi)  F  +  a  fi“6 


„  T z+iu>  fit 
fi2  e  °  dfi 


(k4  -  a4fi2)  T  +  a4fi20 


(3.7) 


(3.8) 


To  evaluate  integrals  in  Equations  (3.7)  and  (3.8),  one  can 
extend  the  variable  fi  into  the  complex  plane  and  use  the  method  of 
Cauchy  contour  integration  for  their  evaluation.  The  choice  of  the 
branch  cuts,  contour  integration  paths,  etc.,  are  enumerated  below. 


3.1.1  Choice  of  branch  Cuts.  The  existence  of  the  function 


2  2  2  2 

F  =  A  -  a  fi  /h  in  the  integrals  requires  that  integrand  be  made 
single-valued  functions  before  Cauchy's  theorem  is  to  be  applied.  This 
requirement  can  be  satisfied  by  introducing  branch  cuts  in  the  complex 
plane.  The  branch  points  are  defined  by  T  =  0  ,  which  yields  fi)}  =  + 
located  on  the  real  axis.  To  allow  an  effective  application  of  the 
methods  of  operational  analysis,  a  complex  value  of  the  wave  number 
k  =  k^  +  ic  (where  C  <<  1  )  has  to  he  used  instead  of  real  valued  k 
in  order  to  avoid  difficulties  in  constructing  the  paths  of  integration. 
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?  ?  i/2  _  1/2  _  1/2 

r  =  (k“  -  Cl  :r)  =  (k  -  a ft)  (k  +  aft) 


where  each  terra  associated  with  one  of  the  pair  of  branch  points 
k  — 

ftQ  =  +  -  .  Substituting  (k  -  aft)  =  r^expCiO^  +  iu)  and 
(k  +  aft)  =  roexp(i02)  into  f  ,  one  obtains: 


_  1/2  _  1/2 
r  =  (k  -  aft)  (k  +  aft) 


1/2  ri(0  +  0  +  TK 

(rir2)  exp  -  2 - 


=  r  exp(i0  ) 
o  o 


where 


ro  *  ‘'iV 


and  0 


ex  +  e2  +  tt 


are  the  amplitude  and  phase  angle  of  the  polar  representation  of  T  , 

respectively.  Thus,  the  multivalued  angle  function  0  can  be  made 

o 

single-valued  on  the  two  Riemann  sheets  defined  by: 


top  Riemann  sheet 


0  +  6  +  IT 

1  2 


bottom  Riemann  sheet 


61  +  62  -  ” 


where  the  positive  direction  of  0^  and  0^  is  measured  from  positive 
real  axis  of  the  complex  plane  ft  in  the  counterclockwise  sense  as 
shown  in  Figure  3.2. 


3.1.2  Evaluation  of  Integrals  and  •  Examining  the 

integrands  in  and  P2  together  with  outward  radiation  conditions 

ImT  >  0  and  ImT  <  0  ,  respectively,  it  can  be  shown  that  the 
integrations  must  be  performed  on  top  Riemann  sheet.  The  closed  paths 
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of  integration  must  be  in  the  fourth  and  the  first  quadrants  of  the 
complex  fi-plane  for  P  and  P.,  ,  respectively,  as  shown  in  Figures  3.3a 
and  3.3b. 

The  closed  path  for  the  integral  may  be  written  as: 


_  I’z-iio  Qt 
e  ° 


dfi 


(k  -  a  Si  )  F  +  a  Si  p 


=  — 27T i  {  I  Residues  } 


o  c  — 
R1 


-2iri{  E  Residue  - 


cR^  -i°° 


9  Fz-iu>  fit  1 

Si  e _ °  dSl _ 

—A  42—  42 

(k  -  A  )  r  +  a  sr$ 


(3.11) 


Similarly,  the  integral  P,  may  be  written  as: 


r ,  =  +2ni{  E  Residues  1  ^ 


-  r  z+iu)  si t 

SI2  e 


d  SI 


c  i°° 
R2 


(k4  -  a!* Si")  F  +  cAt3J 


(3.12) 


3.1.3  Evaluation  of  the  Pole  Contributions .  The  poles  of 
integrals  P ^  and  P0  are  defined  as  the  homogeneous  roots  of  their 
denominators : 

( k '  -  a* Si2)  T  +  cM220  =  0  .  (3.13) 

1/2 

—  2  — 9 

Substituting  otfl  =  (k~  -  F  )  into  Equation  (3.13),  one  obtains  the 
following  cubic  equation: 

—3  —2  T  --4  2—2  —2 

F  -  3r“  +  (k  -  a k  )  +  k“3  = 

a 


0 


(3.14) 
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-  _2  2  ,  1/2 

The  multivaluedness  of  T  =  (k~  -  a  £2")  indicates  that  there  are  six 

complex  Q-poles.  Numerical  computation  of  these  roots  shows  that  four 
of  those  complex  poles  fall  on  top  Riemann  sheet  (Imr  >  0)  and  two  real 
poles  fall  on  bottom  Riemann  sheet  (ImF  <  0).  However,  only  the  two 
poles  on  top  Riemann  sheet  that  fall  inside  the  closed  contour  of 
integration  are  needed  in  the  Residue  theorem.  Let  the  poles  that  fall 
inside  the  contour  of  integration  of  P  and  be  designed  as 

ftp^(ReT  >  0,  ImT  >  0)  and  f2^n(Rer  >  0,  ImF  <  0)  ,  respec  tively . 

One  can  approximate  the  position  of  the  pole  that  falls  inside 
the  paths  of  integration.  The  approximate  expressions  for  Sip  and  Tp 
are : 

—2 

ftp  %  ~  ,  (3.15a) 

a" 

which  represents  the  frequency  spectrum  for  an  unloaded  Euler  plate 
model,  and 

—  'u  k*  2  a, 

Fp  ^  +i(-2-k)  ^  +  i  —  ,  (3.15b) 

a 

respectively.  These  approximations  are  good  for  large  values  of  k  . 
Closer  examination  of  these  expressions  shows  that  they  are  valid  for 
frequencies  above  the  coincidence  frequency  where  the  acoustic  mass 
loading  is  minimal.  Comparing  these  approximate  values  with  exact 
values  computed  from  Equations  (3.13)  and  (3.14),  in  Figures  3.4  and  3.5, 
one  can  conclude  that  these  are  very  good  approximations  for  k  >  1.0  . 

3.1.4  Computation  of  Residues.  The  residue  of  the  integrands 
at  the  simple  poles  are,  respectively: 
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o 


Substituting  for  I  (k.z.t)  from  Equation  (3.20),  one  obtains: 


iSu^e 


I(r,z,t) 


2  Th 


{-I1(r,z,t)  +  I2(r,z,t)}  H(cQt  -  z)  ,(3.22) 


where 


r^r.z.t) 


3a6n“  -  2a4k2  -  a4!4  +  2'T,gf 

P1  pl 


(3.23) 


and 


I?(r,z,t)  = 


,k  exp (T  z  +  iu  o  t}  J  (kr)dk 
"  o 


"p2  p2"  r>2  o"p2 _ 

3a6D20  -  2a4k2  -  +  2a43f 


P-^ 


(3.24) 


It  should  be  noted  that  the  degree  of  the  polynomials  in  the  numerator 
and  denominator  of  the  integrands  of  Equations  (3.23)  and  (3.24)  are 
equal.  The  two  must  he  divi.  d  so  that  the  remainder  function  has  a 
denominator  that  is  at  least  one  degree  higher  than  the  numerator  and 
consequently,  Hunke  inverse  transformation  would  converge.  Dividing 
out,  one  obtains: 


5  „  7T 


i  (r,z,t)  =  N  *  +  -1 — IjpJL' ilL_L121_jel 

1  1  l  .  >  9  ,5/-.  0-2  _  4,-2  2-4  4  - 

-■  ‘  (  1  ,  ~  2  a  k  -  ak  +  2a  ST  ) 

P1  pl 


X  ex^rpl2  -  i"onplt)  kJQ(kr)dk 


(3.25a) 


and 


I_(r,z,t)  = 


f  /-i  ^  2ct\2rp2+i[3a6-p2-2a4k2-a2^+2aAgrP2] 


f-i-  +  _ HLil2_  _ 

4a5  2a5(3aV9  -  2a V  -  a2^4  +  2a4 

n  A 


6]p2> 


X  exp(r  z  +  itJ  n  „t}  kJ  (kr)dk 
p2  o  p2  o 


(3.25b) 
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l^r.z.t) 


I  (r.z.t)  +  T  (r.z.t) 


(3.26) 


10  (r.z.t)  =  I  ,  j  (r.z.t)  +  I.,.,  (r.z.t) 


(3.27) 


Iu(r.z.c) 


!  2  ^  r  * 2  ’ c ) 


— “T  k  J  (kr)  exp{T  z  -  ico  SI  ,  t}  dk  ,  (3.28) 

.  j  o  pi  o  pi 

-a  J 

oj  o 

r  2aV.  T  -  i[3a6S72  ,  -  2a4k2  -  a2k4  +  2a4gf  ,] 
J _ Pi-Jll _ Pi _ _ Pii 

I1  2a  (3a  ft.  -  2a  k  -  a  k  +  2a  3r  ,) 

n  P  1  P1 


k  J  (kr)  exp{F  z  -  iw  f2  t)  dk  ,  (3.29) 

o  pl  o  pi 


1,  .(r.z.t)  =  — r  Ik  J  (kr)  exp[T  z  +  iui  SI  t  ]  V  dk  (3.30) 

-1  l  °  P 2  °  P2  J 


I  (r.z.t)  = 


f2a5Sl  .,r  „  +  i [ 3a6Q2.,  -  2a4k2  -  a2k4  +  2a4SrP2^ 
J  .  .  P~  P~ _ P~  _ 1 _ 1 

2a5(3a6Q2  -  2a4k2  -  aV4  +  2a4t3f  .,)  > 

p2  p2' 


*  k  J  (kr)  exp{F  z+  ito  f)  t)  dk 
o  p  2  op.: 


(3.31) 


3.2.1  Evaluation  of  Integrals _ 1  ,  and  _l. In  the  evaluation 

of  the  integral  in  Equations  (3.2l>)  and  (3.31),  it  is  advantageous  to 
change  from  Bessel  functions  to  Hankel  functions  in  the  integrands. 

Thus,  use  is  made  of  the  identity: 

J  (kr)  =  A-  [ll(l)  (kr)  +  H(;)  (kr)  ] 

O  l  O  o 


H(1) (kr) 

o 


28 


Lott ing 


. ( i) » (2) .  *  AT 


(kr)  ^  77kr  exp{+i(kr  -  £)} 


which  is  the  asymptotic  form  of  the  llankel  functions  for  large  arguments 
(kr  >  1)  ,  aiul  neglecting  higher  order  terms.  Equations  (3.29)  and  (3.31) 
can  he  written  as: 


I12(r,z,t)  j  2iir 


—  f  [2ct5S2pirpl  -  i(3ct6fip1  -  2q4k2  -  qV*  +  2q4grp]  1 


2u5(3a6fi2  -  2a4k2  -  aV*  +  2a4Br  ,) 

pi  pi 


x  k^“  exp(F  .z  -  iw  S2  ,t  +  ikr  -  7-}  dk 
pi  o  pi  4 


(3.32) 


j—  {2a5f2  T  +  i[3a6fl2_  -  2a4k2  -a2k4+  2a43T  0ll 

I  (r  z  t)  =  /-i-  - - p2 _ E2_ 

■)v  ’  ’  ;  /  2nr  5  6  9  4- -?  2—4  4  — 

2a3(3aV.,  -  2a  V  -  aV  +  2a  gr  ) 

p2  p2 


x  kk/,“  exp{T  0z  +  iai  fi  t  -  ikr  +  i-r-}  dk  .  (3.33) 
p2  o  p2  4 


Equations  (3.32)  and  (3.33)  are  written  in  a  form  that  permits 
evaluation  by  the  saddle-point  method.  Tt  will  be  found  convenient  to 
transform  the  observer  coordinates  l rcm  cylindrical  (r,z)  to  spherical 
(R,iJ>)  coordinates  by  the  relations  z  =  Rcosiji  and  r  =  RsiniJ)  ,  where  R 
is  the  radial  distance  from  the  origin  and  i(i  is  the  polar  angle 
measured  from  the  z-axis.  Equation  (3.32)  can  be  written  in  the  form: 


nr  oo 


i12(R,i|>,t)  ^  27TKsini|) 


4  1  5 
e  n 


R*t>  (k) 

FT(k)e  1  dk  ,  (3.34) 
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where 


Ft(k)  = 


k1/2{2a5Q  T  .  -  ipaV.  -  2a4k2  -  a2k4  +  2aV  .J} 

_ pl  pl _ pl _ _ _  pi 

(3a6rr.  -  2a4k2  -  a2^4  +  2a4  gf  .) 

Pl  pl 


and  the  term  in  the  exponential  function  is  defined  by: 

10  t 

$T(k)  =  T  cos4>  -  i— §~0  +  iksin$ 

L  p  1  K  p  1 

The  saddle  point (s)  k^  of  the  above  integral  is  (are)  defined  as  the 
value(s)  of  k  for  which  the  first  derivative  with  respect  to  k  of  the 
exponential  term  ^^(k)  is  zero,  i.e.. 


d$(k) 

dk. 


0)  t 
o  . 

R  '"pl 


=  (P  ,  COS'P  -  i-§— P  ,  +  1  sinij) 


=  0  . 


(3.35a) 


k=k 


If  one  substitutes  Equation  (3.15)  for  and  F^  into  Equation  (3.35a), 

one  obtains  a  fourth-degree  polynomial  equation  whose  roots  are  the 


saddle  points  kg  as  follows: 

?  2  9  o 

r4,Acos6  AT  ,  .  —  3,Asin$i,  .  —2,  Acos  ij>  sin  <t>  ,  AT  , 

k  ( — ^ - +  k  1 — ~J~>  +  k  1 - - T*  +  — 2~' 

a  a  a  a  a  a 

+  k  (-  4T~"^)  +  1.0  =  0 


(3.35b) 


The  saddle  points'  locations  in  the  complex  plane  are  functions  of  R/h  , 

$  and  T  .  Thus,  as  T  or  R/h  decrease,  the  saddle  point  moves  far 
from  the  origin.  It  can  be  shown  that  only  three  of  the  roots  of 
Equation  (3.35b)  are  saddle  points. 

The  path  of  steepest  descent  (PSD)  is  defined  as  (a)  ImF^(k)  = 
Imr^.(ks)  ,  and  (b)  the  value  of  the  exponential  is  largest  at  the  saddle 
point  and  decreases  most  rapidly  along  the  path.  From  numerical  analysis, 
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the  proper  path  of  steepest  descent  (PSD)  and  saddle  point  are  shown  in 
Figure  3.6. 


The  integral  of  Equation  (3.34)  can  be  written  in  the  form: 


R$T (k) 

FT(k)e  dk  + 


R$r(k) 

F^(k)e  dk  =  2iri  I  Residues 


PSD . 


Since  no  poles  are  located  inside  the  closed  path  of  integration, 
therefore,  E  Residues  =  0  ,  one  obtains: 


R1>(k) 

Fj(k)e  dk  = 


R4>(k) 

Fj (k) e  L  dk 


PSD 


To  evaluate  the  integral  in  the  PSD,  either  one  of  two  methods  must  be 
employed.  Since  the  integrand  of  Equation  (3.34)  has  poles  located  at 
k  ($)  ,  the  integrand  may  become  very  large  when  the  saddle  point 
approaches  the  pole.  As  the  distance  R/h  or  the  time  T  increases, 
the  saddle  point  approaches  the  pole.  Hence,  the  regular  saddle  point 
(RSP)  method  would  be  appropriate  if  the  saddle  point  is  located  far  from 
the  pole.  If  the  saddle  point  is  located  near  the  pole,  then  the  modi¬ 
fied  saddle  point  (MSP)  method  is  appropriate.  Thus,  for  such  distances 
where  the  saddle  point  is  far  from  the  pole,  i.e.,  T  ^  cosij)  and 
R/h  <  100  ,  the  solution  can  be  obtained  as  follows: 


/ — 3^ -  .5  R-fT(k  .) 

=  /— — — - — —  J2 —  F  (k  )e  ^ 

/  n  ,fi"  n.  \  Z  r  t  v  1  > 


R*r(k) 

FJ.(k)e  dk  - /R(j>'.(k  )  5  ‘Iv"sl 

v  1  s 1  2a 


PSD 

and  hence, 


(3.36) 


/ — 1 - h5  ~i7r 

I12(R^,t)  =y  27iRsin<^  ^5  °  4 


/  -2 


Frr(k  ,)e 


R^(ksl) 


R4>"(k  .)  IT  si' 
I  si 


valid  for  0  >  0  . 


(3.37) 


Path  of  Integration  anil  Path  of  Steepest 
Descent  (PSD)  in  the  k-Plane 
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Following  the  satin.'  procedure  for  the  evaluation  of  the  integral  1  , 

one  obtains: 


i  Ti 
5  A 


I0  ,(K,4>,t) 


/  I  hV  f 


'  O  (\  \ 

'  2» K“  •»*  2»V  1  ‘  /“+i'|  <Vb2>  "  52 


Mn(k»d 


(3.3H) 


When  the  saddle  point  k^  comes  close  to  the  pole  k^  ,  the  saddle-point 

method  must  be  modified  to  take  into  account  the  contribution  of  this 

pole.  Jones  [21)  has  outlined  a  procedure  for  evaluating  this  tvpe  of  an 

integral  when  the  saddle  point  is  near  a  pole.  The  procedure  is  to 

isolate  the  pole  k  that  occurs  near  the  saddle  point  ,  and  then  to 
P 

approximate  the  remainder  bv  its  value  at  the  saddle  point  k 

s 

Following  Jones'  procedure  for  the  evaluation  of  the  integrals  T j and 
l ,  ,  one  obi  a  ins  : 


and 


-  i  'I 
A 


,  2Ksin,s  ;  T  |  •ll'l(ks)eXplR't’l(ks)  +  2  S  1 


x  orfe{-(  K)  '  ^  7.^  } 


I-0(K,  ,t) 


,3  T 

It  A 


2Rsin<)>  5 
/  2n 


(3.  Id) 


1  ’  1 

ihi  !  (k  ;)exp|  lM’n  (k.)  + 
erfc(-(V/-7.11l  ,  (3.  AO) 


whore 


ht(k  )  (k  -  k  )  H, 
F  s  s  |-  1 


h,  r(k  )  =*  (k  -  k  )  V 


1  1 


P  11 


1/2 


7,  -  (k  -  k  )  [ <t»" ( k  )  ! 

I  p  s  I  s 


Z1T  “  (k  -  k  (A  (k  3 
IT  l>  s  1 1  s 


Vk  ) 

l  s 


,co«(|>  -  i  l  ft  -4  ik  slue}' 
pi  o  pi  s 


W 


1'  cosiji  +  ii  i:  ,  -  ik  sinA 
1’2  o  p.’  s 


orlc(x)  tin-  comp  1  i'iikmi  t  .u  v  error  I  unction. 


o  K 


the  saddle  point  »  k 


^’(kj 


il  (k) 


d‘  l', ,  (k) 


lor  l.trjte  values  ot  l  hi'  argument  i'l  the  error  funot  ion.  the  simplit  ied 
expressions  in  r<|uations  (1.173  and  (1.1S3  can  ho  used  instead  ot 
Equations  (3.  ll>)  and  (3.40). 


3.2.2  Evaluation  ol  liitoprals  l  and  1.  .  Using  tlradshtevn 

and  Ryrliik  [32]  table  of  integrals  to  evaluate  the  integrals  1  ;  and 
T  ,  |  .  '"'il  left  inp,  ::  ■  RcosiJ)  and  r  RsiniJ>  in  Equal  ions  t  1.383  and 
(  1 .  10) ,  one  oht  a i ns  : 


1 1  I  (r  ,ij>,  t  ) 


,  kd  (kRsinl')  exptl'  .Rcesi}*  -  in'  s'  ,tl  dk  . 
'  "  P'  >’  I'l 


Rubst  itut  inp,  tor  S'  and  1  hv  their  approximate  values  and 

pi  p !  11 

eva  Inal  inp,  the  integral,  one  oht  a  ins: 


I  (R,4>,t) 


4  /  r  ,  2 

Lh  I  a _  11  gRsln  <J> 

5  \im  -  cos$]  *  x*  2  All  l  r  -  cos<|>] 


.  <t>  >  O. 


0.41) 


l21(K,s\t) 


k.l  (kUsin>J>)oxp{r  .  Roos^  +  i to  Si  ,1  ]  dk 
5  o  p  2  o  p 7 


-ih  I  g _  t  it  gRain 

2a5  t-Mi  -  cos*!  1  [  3  4h[T  -  cos^] 


<3. A2) 


3.2.3  Final  Solution  oi  t  ho  ln(  egra  1  1  (R  ,  1 1  .  Subslitut  ini; 

lunations  (307),  (3.3S),  (3.40,  and  (3.42)  into  Filiations  (1.2t>)  and 
(1.27),  and  substituting  t  ho  resulting  equations  into  Fquat ion  (3.72), 
t  lion  the  f  inal  simplified  solution  becomes: 


I  (K.sOt) 


„  .  •  ,  i 

o  .  'iksin  ■;  I  \ 

5  sin  7.  ..  |T  R(T  -  cos<p)  4* 

.  . ,  'di  i  I  -  ooso  i  , 

4nRh  (l  -  cos^) 


IOC  (  ■< 

+  °  )  ■'  ---  F  (k  1  exp  |  R'J'.  (k  )|  + 

4-iRh.sin7  1  ,) 


F  (kj  exp[!On(ks)  I  ■  lt(T  -  cos*) 


(3.40 


Similarly,  substituting  Fquat ions  (1.30),  (1.4l1),  (  1 .  1  O  ,  and  (1.42) 
into  liquations  (3.26)  and  (3.27),  and  suhst  i  t  ut  i  ng  the  result  (hr 
equations  into  Fquat  ion  (1.27),  then  the  tinal  solution  becomes : 


■  ..  7«7, 
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and 

0  =  arc  tan  ($"/$") 

3.2.4  Special  Case  and  Approximate  Method .  In  the  case  when  the 

observation  point  is  on  vertical  axis  <j>  =  0  ,  the  Hankel  function 

(1)  (2) 

H  “  (kRsin<J>)  would  become  infinite.  To  avoid  this  problem,  one  can 
substitute  J  (kRsinif)  =  1.0  in  Equations  (3.29)  and  (3.31),  and 
evaluate  the  radiated  pressure.  For  detailed  solution,  see  Appendix  A. 
For  the  approximate  solution  of  the  radiated  pressure  from  an  infinite 
elastic  plate,  see  Appendix  li. 


CHAPTER  IV 

NUMERICAL  RESULTS 

4.1  Introduction  and  Input  Data 

Numerical  evaluation  of  the  radiated  pressure  at  different 
observation  angles  i}>  and  nondimensional  distance  R/h  and  time  T  is 
performed.  Since  the  transient  acoustic  pressure  is  a  function  of  all 
of  these  variables,  one  can  only  plot  the  pressure  at  a  specific  value 
of  R/h  and  $  as  a  function  of  time  T  .  To  investigate  the  influence 
of  one  of  these  parameters  on  the  total  time  signature,  the  other 
parameters  are  kept  constant.  The  solution  for  the  acoustic  radiated 
pressure  time  signature  for  different  load  histories  are  also  presented 
in  this  chapter.  In  all  of  these  plots,  the  variation  of  the  pressure 
as  a  function  of  time  is  very  rapid  for  times  just  after  the  first 
arrival.  This  made  it  very  difficult  for  the  computer  plotter  to  plot 
a  smooth  sinusoidal  curve  as  predicted  by  the  solutions.  To  avoid  this 
problem,  the  relative  maxima  and  minima  were  connected  by  straight  lines. 
The  resulting  curves  look  smooth  near  the  first  arrival,  hut  become 
saw-tooth  shaped  for  longer  times.  However,  the  location  of  these 
relative  maxima  and  minima  are  known  from  the  solutions,  such  as  those 
given  in  Equation  (3.44).  Thus,  the  connected  points  fall  on  the 
envelope  of  the  smoothly  predicted  curves. 

For  the  sake  of  numerical  computation,  the  plate  was  taken  as  a 
steel  plate  immersed  in  water.  The  relevant  physical  data  are: 


■ 
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E  39  *  10*1  psi,  or  20.68  *  10^  N/nT'  , 

p  -  7.28  *  10*  — 5 ~r~~  .  or  7.6S  kg/m'  , 

in 

■> 

p  =  0.94S  x  10*  — ,  or  1  kg/m'  , 

o  .  4 

in 

c  =  6.0  *  10*  in/sec  ,  or  1500  m/sec 

o 

and 

v  =  0.3 

4.2  Radial o-l  Pressure  for  l)i  fforonr_  Typos  of  Appl  ied  Foroo 

4.2.1  Impulse .  Numerical  values  for  the  radiated  pressure  of  an 
infinite  elastic  plate  excited  by  an  impulsive  point  force,  for  different 
observation  angles  and  distances  were  obtained  by  using  the  above  data 
and  'liquations  (3.45),  (3.46),  and  (A.  4). 

In  order  to  examine  the  relative  cent  filiation  of  the  first  and 
second  terms  in  liquations  (3.45)  and  (3.46),  their  numerical  values  wore 
plotted.  Thus,  plots  of  the  contribution  of  the  first  term  for 

O  P 

R/h  =  50  and  $  =  30  and  i}1  =  60  are  shown  in  Figures  4.1  and  4.2, 
respectively.  Plots  of  the  numerical  contribution  of  the  second  term  of 
these  equations  are  similarly  plotted  in  Figures  4.3  and  4.4  and  the  sum 
of  the  two  terms  shown  in  Figures  4.5  and  4.6.  In  order  to  examine  the 
influence  of  K/h  ,  plots  of  the  impulse  response  are  shown  for  R/h  =  100 
and  1 000  in  Figures  4.7  and  4.8,  respectively.  To  examine  the  oscillatory 
nature  of  the  impulse  response  for  long  times,  the  time  scale  was  expanded 
am!  t lie  result  plotted  in  Figures  4.°  through  4.11  for  v*1  =  30  and 
R/h  =  50,  100  and  1000,  respectively.  Using  Equations  (11.8)  and  (B.9) 
for  the  solution  (RSP)  ,  Equations  (11.13)  and  (11.14)  for  the  solutions  MSP, 
numerical  results  wore  obtained  for  radiated  pressure  predicted  from  those 


Magnitude  of  the  First  Term  of  R) 


eco 


P(2~h  R) 


P(2rh  R) 


1.1  2.3  3.3  4.7  3.9 


T 


l’ip.ure  4.10  Plot  of  tlu’  Acoustic  Unit  Impulse  Response  with  an 
Expanded  Time  Scale  for  the  Observation  Point  at 
4>  =  30°  and  R/h  =  100 
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approximate  solutions.  Figures  4.12  through  4.15  and  Figures  4.16 
through  4. Id  show  the  impulse  response  of  the  radiated  pressure  versus 
time  for  ditferent  observation  angles  and  observer  distances  for  MSP  and 
RSP  solutions,  respectively. 


4.2.2  Applied  Force  with  a  Const  ant  _Mar.n  it  ude  .  Such  an  applied 
force  is  defined  by  Q(t)  =  f^H(t)  .  Substituting  applied  force  and 
Equation  (3.43)  into  the  convolution  integral  in  Equation  (3.3),  and 
neglecting  the  very  small  contribution  of  the  second  term  of  Equation 
(3.43),  one  obtains: 


when 


and 


t  N ,  s  i  n 
f  1 


P  ( R ,  O ,  t )  =  f 


arc 


4"Rh‘ 


otRs  i  n  <P 
4h 


N-2 
i  c 


T  - 


-  cosO1 


re 


rc 


T  - 


o  1 

-  cos>?j 


H(T 


-  cos<fOdT 


(4.1) 


1’ransl  erring  the  integration  variable  t  to  u  bv  lettii 


N2 

Tc 

o 

E 


and  du 


c  N ,  dT 

o  2 _ 

1  c 

R(F  -  -  cos*)' 


Equation  (4.1)  can  be- written  as: 

N, 


F  ( R ,  s' .  t ) 


v  i>  T-eos(J> 

- 1  NR  r 

o  1  I  sin  u 


H(T  -  cos*,*')  du  . 


(4.2) 
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Figure  4.12  Plot  of  the  Aeoustie  Unit  Impulse  Response  with  an 
Expanded  Time  Seale  for  the  Observation  Point  at 
<J>  =  30°  and  K/ It  =  30  Obtained  by  the  Modified 
Saddle-Point  Method  (MSP) 
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Plot  of  the  Acoustic  Unit  Impulse  Response  with  an 
Expanded  Time  Scale  for  the  Observation  Point  at 
<}>  =  60  and  R/h  =  100  Obtained  by  the  Modified 
Saddle-Point  Method  (MSP) 
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To  evaluate  the  Integral  of  Equation  (4,2),  one  can  substitute 
sin  u/u  ,  with  its  power  series  expansion,  and  then  integrate  term  by 
term.  Letting 


-  I 

„  ill  '  1 


+  D! 


and  substituting  into  Equation  (4.2),  and  evaluating  the  resulting 
integral,  one  obtains: 


P(R,4>,t)  = 


f  a  3  00 

l 

4iTh"  n=0 


-1)  = - = - r 

r  -  cosi>) _ 

(2n  +  1  )\ ,!n  Vl)  ! 


•  (4.3) 


The  normalized  radiated  pressure  from  an  infinite  elasLic  plate  excited 

by  a  point  force  Q(t)  =  fll(t)  is  then  given  bv: 

o 


2  V  i 

n=0 


cu-r^i 

-  1  -  cosO  - 


l)(2n 


-  cos  ) 


(4.4) 


Figure  (4.20)  and  (4.21)  are  the  plots  of  radiated  pressure  versus  time 

for  applied  force  Q(t)  =  f  H ( t )  ,  R/h  =  50°  and  it>  =  30°  .  The 

o 

normalization  pressure  P  is  the  pressure  computed  for  an  infinite 
elastic  plate  excited  bv  a  time  harmonic  point  force  of  amplitude  f 

o 

at  coincidence  frequency  and  evaluated  at  (lie  observer  angle  4*  =  0° 
given  by: 

f  otfj 

|P  |  = 

o 


2n(R/h)ir 
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'  App  1  i e d  Square  Pn  1  sc  Force.  Such  an  applied  force  can  be 
represented  by  Q(t)  =  f^Kt^  -  l)  .  Substituting  applied  force  and 
Equation  (1.43)  into  the  convolution  integral  Equation  (3.3),  and 
neglecting  tin'  very  small  contribution  of  the  second  term  of  Equation 
(3.43),  one  obtains: 


P(R,*,t)  =  f  N, 


|T - £ - cost}) 


[T - ^ - cost,1)] 


H(T - - - cos<|>)  dx  , 


for  T  <  T  ,  (4.5) 


P(R,v\t)  =  f  N  I  - 

o  .1 


[T - J.-  ~  cosj 


- — -  H(T  -  —  -  cos*)  dx  , 

tT  -  ~  -  cos*] 

for  T  >  T  ,  (4.6) 

o 


where  T  = 
o 


The  Integral  in  Equation  (4.5)  is  similar  to  the  integral  in 
Equation  (4.1),  therefore,  following  the  same  procedure  t:o  evaluate  the 
integration,  one  obtains: 


abf  00  (-1)" - - - 1 

m  H« . •»  • 


for  T  <  T  .  (4.7) 

—  o 


E°r  T  <  T  ,  the  normalized  radiated  pressure  from  an  infinite  elastic 

plate  excited  by  a  point  force  Q(t)  =  f  ll(t  -  t)  is  then  given  bv: 

00  • 


n  r  No  1 2n+i 

a,  (_nn  2 

_L_  _  V  -T  -  cos*J 

Po  ~  2  h  n:0  2("  +  l)(2n  +  iTT-*1^  '  cos<|,)  'forTi  T0 


(4.8) 


Transferring  the  integration  variable  T  to  u  by  letting 


and  du  = 


T  -  —  -  cos  * 


c  N  dT 
o  2 


R(T - -  -  cos*)‘ 


Equation  (4.6)  can  be  written  as: 


T-T  -cos* 

p(R,^,t)  =  foNl  ° 

n2  ° 

T-cos* 


T  >  T  , 
o 


(4.9) 


P(M,t)  = 


f  NR 
o  1 


T-cos* 


sin  »  „/  2. 

~  H(ir)  du 


f  N  R  T-T  -cos*  .  N 

i  °  1  0  sin  u  2,  . 

+  co  u  "<T>  du  •  T  "  To  • 


(4.10) 


Hie  integrals  of  Equation  (4.10)  are  similar  to  the  integral  of 
Equation  (4.1),  therefore,  following  the  same  procedure  to  evaluate  the 
integration,  one  obtains: 


471b2  'n-0  (2ll+  l)(2l'  +  l)!  [l* 

?n+l 

No] 

x  H  (T  -  cos4>)  -  - = — ‘ - r  H(T  -  T  -  cos  A) 

1-1  -  coso  o 

°  J  J 

for  T  >  T  .  (4.11) 

o 

The  normalized  radiated  pressure  from  an  infinite  elastic  plate  excited 

by  a  point  force  t)(t)  =  f  H(t  -  t)  is  given  bv 

o  o 


(  00 

„  11 

N 

y 

(-1) 

n2 

i  *• 

K  n=0 

(2n  i  1 l (2n  +  1) ! 

T  -  co so 

K 

2n+l 

1 

- = - 

H(T  -  T  -  cost})) 

T  -  T  -  cosO 

o 

o 

j 

for  T  N  T  .  (4. 12) 

o 

1  iguros  4.22  and  a. 23  are  the  plots  of  radiated  pressure  versus  time  for 

applied  force  Q(t)  =  f  lift  -  t)  and  T  >  T  for  T  =  cost  -  0.1  an 

o  o  o  o 

"  0.5  cos ^  ,  respectively,  the  forcing  frequency  =  1.5  , 

H/b  ---  SO  and  4>  30°  . 

4.2.4  Applied  CU'-l'ulse  Force.  Such  an  applied  force  can  be 
represented  by  Q(t)  =  f^cosijjtllCt  -  t)  .  Substituting  applied  force 
and  liquation  (3.43)  into  the  convolution  integral  liquation  (3.3),  and 
neglecting  the  very  small  contribution  of  the  second  term  of  liquation 
(3.43).  one  obtains: 
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x  H(T  -  T  -  cost}))  - 


cos  (oj  a)  cos  (id  — )  -sln  — 

1  1  U  u 


T-cosij) 


iin(coia)sin(wi  — -)  — ^  u  duj  H(T  -  cost?) 


for  T  >  T  ,  (4.15) 


where  a  =  T  -  cos<i>  . 


Equation  (4.15)  can  he 


;  rewritten  as  follows: 


P(R,.r\t)  -  |  [P  +  P  +  P  _  Pi  j  h(T  -  cos«) 


CP31  +  P32  +  P41  “  Pa-]  H<T  -  T„  “  cos«fr)f,(4.16) 


cost)  a  T-cos4>  n 

— r~  s  i  n  ( u  +  — 3.)  — . 

u  •  U 
J 
o 

N2 

costa,  a  [  T-cos*  N 

_ _ L_  •  /  3.  du 

'  sm(u  -  — )  — 

C  U  U 


sinw  a  T-cosij)  N 

_ 1  ,  3  s  du 

2  cos (u  - )  — 

L  u  u 


(4.17) 


(4. IS) 


(4.19) 


smw  a  T-cosif 


/  .  3.  du 

cos  (u  +  — )  - 

u  U 


(4.20) 


o 


I 


—  1-1  -cossp  .. 

coso),  a  o  N  . 

1  .  ,  ,  lx  du 

- ^ —  sin(u  H - ;  — 

2  u  u 


(A. 21) 


cosco  a  1-1  -cos4>  N  . 

1  o  .  ,  3,  du 

- ^ —  sin  (u  -  — )  — 


(A. 22) 


sinu),  a  1-1  -cosij)  N  , 

1  o  t  3,.  du 

- - -  cos  (u - )  — 

2  u  u 


(4.23) 


sinOJ,a  T-T  -cos^  N.,  . 

1  o  /  ,  3,  du 

-  cos(u  1 - )  — • 

2  u  u 


(4.24) 


1/2  du  . 


LoCCinc  y  =  u  +  —  ,  and  dy  =  (y  -  4N-)  —  into  Equations  (4.17) 

"  u  J  u 


and  (4.20),  one  obtains: 


coswa 

1  sin  v  . 

—  -  ->iy 

l  y  y  -  4N 
1 


(4.25) 


si .nw  a 


/  y2  -  41. 


(4.26) 


where 


N“  +  N.;(T  -  cos40 2 
N2(T  -  cosifi) 


ry  1  /  ^  1 

Lot  tine  v  =  u - ,  and  dy  =  (y  -f  4N.J  "  -  into  Equations  (4.18) 

u  J  u 


ind  (4.19),  one  obtains: 
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cosiii  j  a 


sin  v 


-  /  y  +  an.. 


s  i  iiuij  a 


J  y"  +  4N , 


whore 


N  >  ~  N  (T  -  costy )  ‘ 
N,(T  -  cosifi) 


(4.27) 


(4.28) 


ApPlying  th°  S;1"K'  tyP°  °f  transformations  as  above  to  Equations  (4.21) 
thr°USh  mU'  °bt3lns  ^^rals  identic a  1  to  those  la  Equations 

(4.2s)  tlu-o.iRh  (4.28)  with  different  upper  limits.  Equations  (4.24)  nr 


(4.26)  can  be 


rewritten  as  follows: 
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where  J  (x)  and  Y  (x)  are  the  Bessel  and  Neumann  functions, 
o  o 

respectively.  The  integrals  In  Equations  (A. 29)  and  (A. 30)  can  be 
evaluated  by  numerical  methods.  Equations  (A. 27)  and  (A. 28)  can  be 
r ew r i 1 1  e  n  a  s  foil ows : 

20  J 


cosio^a 


/  2 

^  v  y  +  A  N 
0  3  J 


1  /  y  +  /(N 

200  J 


for  G„  <  0 
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s  into  a 


o v  y  +  /,n3 


COS  V  ^ 
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si mo |  a 
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K  <2 /tO  + 

0  -) 

cos  y 

dy 

5 

v  y  +  AN., 

(A. 31) 


for  c; ,  <  0  . 


(A. 32) 


Since  y  ''  AN  ^  ,  when  y  v>  2(10  ,  the  secoiul  integral  ot  Equation  (*i.  Ml 
can  be  approximated  as  follows: 


sin  v  ,  'O 

d  V  'V 


s  i  n  v  ,  „  2 

*■  dv  ,  I  or  v  ','s  AN , 
v  J  3 


200  /  y  +  AN  200 

which  turns  out  to  be  the  Sine  Integral  function.  The  remaining 


Equations  (4.21)  through  (4.24).  Substituting  these  results  into 
Equation  (4.16),  one  obtains  an  expression  for  the  radiated  pressure 
from  an  infinite  elastic  plate  excited  by  a  point  force 

Q(t)  =  f  cosoi  t  H(t  -  t) 
o  1  o 

The  plots  of  the  normalized  radiated  pressure  versus  time  for 

T  >  T  and  with  T  =0.5  cosd>  and  T  =  cosii  -  0.1  are  shown  in 
o  o  o 

Figures  4.24  and  4.25,  respectively.  The  other  parameters  for  these 
plots  were  chosen  as  £2^  =  U)^/u>0  =  1.5  ,  R/h  =  50  and  ij>  =  30°  . 
An  expanded  time  scale  for  the  results  of  Figure  4.25  is  shown  in 
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CHAPTER  V 


DISCUSSIONS  OF  THE  RESULTS  AND  CONCLUSIONS 

5.1  Discussion  of  the  Impulse  Response  Solutions 

The  numerical  results  for  the  radiated  pressure  signature  due  to 

a  unit  impulse  load  are  discussed  first. 

Equations  (3.45)  and  (3.46)  represent  the  solution  for  the 

radiated  pressure  from  an  infinite  elastic  plate  in  contact  with  an 

infinite  acoustic  medium  occupying  the  space  z  >  0  and  excited  hv 

impulsive  point  force  for  different  observation  angles  and  distances. 

The  expression  for  the  unit  impulse  response  consists  of  two 

terms,  each  exhibiting  a  damped  sinusoidal  oscillation  where  the  period 

varies  with  time.  The  relative  magnitude  of  these  terms  can  be  estimated 

from  the  numerical  evaluation  of  the  contribution  of  the  first  term  as 

shown  in  Figures  4.1  and  4.2  to  the  contribution  of  the  second  term  as 

shown  in  Figures  4.3  and  4.4,  respectively,  where  the  range  R/h  was  50 

and  observer  angles  are  =  30  and  60l  .  These  plots  show  that  the 

ratio  of  peak  value  of  the  first  term  to  that  of  the  second  term  is 
2 

0(10“)  .  Thus,  the  contribution  of  the  second  term  can  be  neglected  and 
solution  is  approximated  by  the  first  term  only  for  all  values  of  T  . 
This  is  evidenced  when  one  compares  the  total  solution  shown  in  Figures 
4.5  and  4.6  with  the  contribution  of  the  first  term  only  in  Figures  4.1 
and  4.2,  respectively.  The  general  character  of  the  solution  is  similar 
when  R/h  =  100  and  1000  and  4>  =  30°  as  can  be  seen  in  Figures  4.7  and 
4.8.  Since  the  largest  contribution  to  the  impulse  response  comes  from 
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the  first  term,  one  can  explore  the  solution  more  fully  hv  examining  the 
first  term  only.  The  first  term  represents  a  damped  sinusoidal  oscilla¬ 
tion,  which  decays  to  2ero  as  T  -*■  00  .  These  oscillations  are  shown 
explicitly  when  one  expands  the  time  scale  as  shown  in  Figures  4.9 
through  4.11  for  4*  =  30°  and  R/h  =  50,  100  and  1000  .  The  period  of 
oscillations  varies  with  the  elapsed  time  T  ,  the  observer  angle  4> 
and  R/h  .  The  period  is  very  short  at  the  first  arrival  of  the  acoustic 
wave  and  becomes  larger  later.  The  time  of  arrival  T  of  the  n1'1  peak 
can  be  deduced  from  the  sinusoidal  term  in  Equations  (3.45)  and  (3.46), 
such  that : 


The 

by 


T 

n 


a  (R/h) isin2(i> 
2  (2n  +  D  'T 


COS'}1 


largest  half-period  between  the 
the  time  difference  between  T 

ti=l 


,  where  n  =  0,1,2, .. . 

last  two  observable  peaks 

and  T  _  as  follows: 
n=2 


.  (5.1) 
is  given 


AT  =  T  -  T 

n=l  n=2 


q(R/h) s in  5 
15n 


The  half-period  AT  is  marked  on  each  of  Figures  4.9  through  4.11. 
Thus,  the  time  interval  between  the  last  two  observable  peaks  of  an 
impulse  response  time  signature  of  an  elastic  plate  may  be  used  to 
obtain  an  estimate  of  the  observer  distance  R  if  the  plate  thickness 
h  is  known. 

In  absolute  time  scale,  AT  in  seconds  is  given  by: 


AT  =  —  AT  =  —  — r- 
c  30nc  h 

o  o 

a,  2 

where  a  't  1.0  and  the  mean  value  of  sin  f  was  taken  to  be  0.5  . 
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The  last  peak  of  the  impulse  response  occurs  at  time  T  ^  ,  after  which 
no  oscillation  occurs  in  the  time  signature  and  it  decays  at  T  ^  . 

The  computation  of  the  time  signatures  for  the  exact  solutions 
has  been  discussed  above.  The  most  noticeable  feature  of  the  exact 
solution  is  the  prediction  of  the  first  arrival  of  the  acoustic  pressure 
wave.  The  exact  solution  predicts  a  nondimensional  arrival  time  of 
T  =  cos^  <  1.0  which,  if  translated  to  absolute  time,  is: 


c  cc  's'»  c  c 
o  o  o 

This  means  that  the  elastic  wave  in  the  plate  would  have  travelled  at  an 
infinite  phase  velocity.  Thus,  the  first  vibration  induced  radiated 
acoustic  pressure  wave  would  have  emanated  from  a  point  on  the  plate  just 
below  the  receiver  point  (R,if>)  and  would  have  travelled  in  the  acoustic 
medium  a  distance  z  only.  Physically,  the  phase  velocity  in  flexure 
cannot  exceed  the  shear  velocity  of  the  elastic  plate  and  thus  one  would 
expect  the  arrival  time  of  the  first  pressure  wave  to  have  occurred  later 
than  T  =  cosi|>  .  The  reason  for  this  underproduction  of  the  first 
arrival  is  due  to  the  Euler  plate  theory  used  as  a  model  for  the  flexural 
vibration  of  the  elastic  plate.  The  Euler  plate  theory  is  an  approximate 
theory  of  antisymmetric  motion  in  a  plate.  This  approximation  assumes 
that  the  normal  to  the  plate's  mid-surface  stays  normal  after  deformation 
and  that  the  flexural  vibration  involves  no  rotation  of  an  element  about 


its  mid-surface.  These  approximations  leading  to  the  Euler  plate  theory 
gives  one  dispersion  curve  $1  vs  k  given  in  Equation  (3.15a).  Thus,  the 


phase  velocity  v  ^  is  given  by  the  secant  v  =  STJ/k 
without  limit  as  |f)j 


which  increases 


V  CO 


Thus,  the  dispersion  curve  predicts  an 
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infinite  phase  velocity  as  ft  -*■  00  .  Since  the  Fourier  transform  in 
Equation  (2.5)  integrates  over  ft  from  to  00  ,  one  would  expect  an 

error  when  ft  is  very  large  or  when  t  is  very  small.  Thus,  the  first 
arrival  of  the  pressure  wave  would  have  travelled  at  infinite  phase 
speeds  when  the  integration  covers  |ft  |  -*■  <»  .  It  should  be  noted  that 
the  Euler  theory  is  valid  for  wavelength-to-thickness  ratio  larger  than 
eight,  i.e.,  low  frequencies.  A  more  accurate  theory  of  flexural 
vibration  is  the  Tiraoshenko-Mindlin  plate  theory.  This  theory  includes 
shear  deformation  and  rotatory  inertia,  so  that  the  two  frequency  spectra 
of  this  theory  accurately  predict  the  first  two  nonsymmetric  spectra  of 
the  exact  plate  theory.  Thus,  the  Timoshenko-Mindlin  plate  theory  would 
accurately  predict  the  flexural  phase  velocity  as  ft  -►  00  to  be  Rayleigh 
wave  speed.  Thus,  the  vibration  in  the  plate  would  travel  by  Reyleigh 
speed  first  before  travelling  in  the  acoustic  medium.  Hence,  the  first 
arrival  of  the  pressure  wave  would  occur  later  than  T  =  cos*}1  but 
before  the  direct  acoustic  path  arrival  of  T  =  1.0  ,  since  the  Rayleigh 
wave  speed  for  metals  is  higher  than  that  of  the  sound  speed  in  water. 

Approximate  expressions  for  the  radiated  pressure  were  also 
obtained  and  given  in  Equations  (B.P)  and  (B.10).  This  approximation 
parallels  the  development  of  similar  solutions  obtained  by  Magrab  and 
Reader  [20]  and  Stuart  [12].  This  approximation  was  obtained  by  the 
regular  saddle-point  method  and  henceforth  will  be  referred  to  as  RSP. 

A  better  approximation  than  the  above,  using  the  modified  saddle-point 
method  was  given  in  Equations  (B . 14 )  and  (B.15).  These  solutions  will 
be  referred  to  as  MSP  in  the  subsequent  discussions.  However,  the  two 
approximations  predict  similar  time  signatures,  but  the  predicted 


. 


magnitudes  of  the  radiated  pressure  were  different.  The  two  approxima¬ 
tions  predict  an  oscillating  impulse  response  for  T  <  1.0  and  a  purely 
exponential  decay  in  the  response  for  T  >  1.0  .  This  is  clearly  in 
contrast  to  the  character  of  the  exact  solution,  which  exhibits  damped 

oscillations  before  and  after  T  =  1.0  and  up  to  time  T  =  T  ,,  ,  which 

n=0 

is  the  time  of  arrival  of  the  last  peak  of  the  oscillating  impulse 

response.  Furthermore,  the  exact  solution  predicts  a  time  decay  1/T 

for  T  >  T  „  ,  which  decays  at  a  slower  rate  than  the  exponentially 
n=0 

decaying  response  predicted  by  the  approximate  solutions.  In  addition, 

the  approximate  solutions  cannot  predict  the  first  arrival  of  the 

acoustic  wave,  unlike  that  predicted  by  the  exact  solution. 

The  solutions  predicted  by  MSP  are  shown  in  Figures  4.12  through 

4.15  for  R/h  =  50  and  100  and  i}>  =  30C  and  60*  .  The  MSP  solutions 

predict  acoustic  pressure  peak  levels  of  the  order  0(10  ~)  smaller 

than  those  predicted  by  the  exact  solutions  for  all  T  as  can  be  seen 

when  one  compares  Figures  4.12  through  4.15  to  Figures  4.9  through  4.11. 

They  also  decay  faster  than  the  levels  predicted  by  the  exact  solutions. 

The  solutions  predicted  by  RSP  are  shown  in  Figures  4.1b  through 

4.19  for  K/h  =  50  and  100  and  ij>  =  30°  and  60°  .  The  RSP  solutions 

_2 

also  predict  acoustic  pressure  peak  levels  of  the  order  0(10  ~)  smaller 
than  those  predicted  by  the  exact  solutions  for  all  T  .  Comparing  solu¬ 
tions  predicted  by  MSP  and  RSP,  one  can  conclude  that  the  two  solutions 
predict  similar  time  signatures,  approximately  the  same  peak  levels  and 
exhibit  an  exponentially  decaying  time  signature  after  T  N 


1.0  . 
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5 . 2  Discussion  of  the  Radiated  Pressure  Due  to  Different  Applied  Loads 

Using  the  unit  impulse  response  given  in  Equations  (3.45)  and 
(3.46),  one  can  obtain  the  solution  for  an  applied  point  force  with  a 
designated  time  signature  by  use  of  the  convolution  theorem.  In  order 
that  one  can  compare  the  relative  magnitudes  of  the  time  signatures  of 
the  radiated  pressure  for  different  loads,  the  numerical  values  were 
normalized  to  the  radiated  pressure  from  an  infinite  elastic  plate 
excited  by  harmonic  point  force  at  the  coincidence  frequency. 

The  normalized  time  signature  from  a  plate  excited  by  a  constant 
load  has  been  computed  and  shown  in  Figures  4.20  and  4.21  for  R/h  =  50 
and  9  =  ^0°  •  The  response  builds  from  a  constant  at  first  arrival  and 
increases  with  T  in  an  oscillatory  manner  as  seen  in  the  expanded  scale 
in  Figure  4.20.  The  oscillations  are  due  to  the  continuous  vibration  of 
the  plate  after  it  had  been  excited  at  T  =  0  by  the  constant  load.  The 
period  of  oscillation  is  zero  at  the  arrival  time  and  increases  as  time 
increases.  The  reason  for  zero  or  very  short  period  at  the  arrival  time 
can  be  traced  to  the  inaccuracy  of  the  Euler  plate  theory  for  high 
frequencies,  which  contribute  mostly  near  the  first  arrival.  After 
T  -  1.0  ,  the  pressure  increases  to  a  maximum  dynamic  pressure  and  then 
decays  monotonically  as  1/T  for  T  »  1  as  seen  in  Figure  4.21.  This 
is  physically  reasonable,  since  as  T  ,  the  influence  of  the  constant 

load  becomes  quasi-static  and,  hence,  no  dynamic  pressure  results. 

The  response  due  to  a  square  pulse  is  shown  in  Figures  4.22  and 

4.23  for  two  pulse  lengths.  In  Figures  4.22  and  4.23,  the  radiated 

pressure  time  signature  due  to  pulse  lengths  T  =0.5  cos^i  and 

o 

T^  =  co.Si})  -  0.1  are  shown,  respectively.  The  response  after  the  first 
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arrival  at  T  =  cos^  and  before  the  arrival  of  the  pulse  representing 

the  time  when  the  load  is  removed,  i.e.,  T  =  cos<J)  +  T  must  necessarily 

o 

resemble  the  response  due  to  a  constant  applied  load  as  shown  in  Figures 
A. 21  and  4.22.  However,  the  response  drops  immediately  by  approximately 
20  dB  after  the  removal  of  the  load.  After  the  removal  of  the  load,  the 
plate  oscillates  freely  and  the  peak  value  in  the  response  corresponds 
to  the  peak  value  in  the  unit  impulse  response  curve.  The  response  for 
long  tines  T  again  decays  as  1/T  for  T  >>  1.0  as  clearly  shown  in 
Figure  4.23. 

The  normalised  radiated  pressure  time  signature  due  to  CW-pulse 

load  for  T  =  0.5  cosd>  and  T  -  cos<5  -  0.1  are  shown  in  Figures  4.24 
o  o 

and  4.25,  respectively.  The  response  after  the  first  arrival  T  =  cosift 
is  a  sinusoid  of  the  same  period  as  the  load.  For  T  >  cos4>  +  T^  , 
which  represents  the  time  when  the  load  is  removed,  the  response  drops 
immediately  as  expected.  After  removal  of  the  load,  the  plate  oscillates 
with  increasing  periods  in  a  manner  similar  to  the  impulse  response.  The 
response  for  long  time  T  again  decays  as  1/T  for  T  >>  1.0  . 

5 . 3  Cone lus ions 

The  prediction  of  the  response  of  an  elastic  plate  to  time 
dependent  point  forces  and  the  resulting  time  signatures  of  the  radiated 
acoustic  pressure  was  the  main  objective  of  this  study.  The  radiated 
pressure  was  predicted  for  a  unit  impulse,  a  constant  load,  a  square 
pulse  and  a  Civ'  pulse  at  different  observer  distances  R/h  and  angles  . 
The  following  conclusions  can  be  drawn  from  the  study  of  the  analytic 
solutions  and  their  numerical  evaluation  as  follows: 


i. 
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a.  The  predicted  time  for  the  first  arrival  T  =  cos<{>  is 
shorter  than  the  acoustic  arrival  T  =  1.0  ,  which  is 
physically  more  realistic  than  the  preceding  studies 
since  the  elastic  waves  travel  in  an  elastic  plate  at 
supersonic  speeds  at  high  frequencies.  However,  a  more 
accurate  prediction  of  the  first  arrival  time  would 

have  been  longer  than  T  =  cosij)  and  shorter  than  T  =  1.0 
if  one  would  have  employed  a  more  accurate  theory  of 
flexural  response  for  the  plate. 

b.  The  time  signature  of  the  radiated  acoustic  pressure 
oscillates  with  an  increasing  period  of  oscillation 
due  to  the  free  vibration  of  the  infinite  plate  after 
the  application  of  the  impulsive  force.  This  was 
attributed  to  the  dispersive  nature  of  the  frequency 
spectrum  of  the  elastic  plate. 

c.  The  time  signature  of  the  radiated  pressure  always 
diminishes  immediately  after  the  removal  of  the  load, 
which  is  physically  reasonable. 

d.  The  time  difference  between  the  last  observable  peaks 
can  be  measured  and  are  directly  related  to  the  observer 
angle  <p  and  the  ratio  R/h  .  Hence,  one  can  get  an 
approximate  value  for  R  from  such  a  measurement  if  the 
average  v.alue  of  the  observer  angle  <}>  is  taken  as  45°. 

5.4  Suggested  Future  Research 

It  was  evident  from  the  analytic  solution  thus  developed  that  a 
Timoshenko-Midlin  flexural  plate  theory  would  have  made  early  time 


pi  i'll  lot  ions  nii'io  .ii-iuir.it  i' .  However,  oven  tills  no  re  .ii-rur.it  i*  thoorv 
won  lit  iii' t  l'i-  v .  1 1  lil  I'l'i.iiisi’  It  do  os  not  ;il  low  ooupl  1  n>*  i>l  tin-  I  I  oxnr.-i  1 
w.ivo  with  lottp,  limit  nal  w.ivi-s.  A  thoorv  tli.it  Includes  such  coup  1  tup, 
would  oltlior  ho  out-  ot  tin*  htp.hor  onloi  .ipproxlniat  o  tlioorlos  ot  pl.ito 
vlbr.itiou  dovolopod  by  Nlmllln  or  t  ho  oxaot  thoorv  ot  elasticity.  to 
oithor  case,  such  coupling  would  accui  .it  high  rroquono  tos ,  which  again 
would  oorrort  tho  short  i  lino  rospouso. 
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APPENDIX  A 


ACOUSTIC  PRESSURE  RESPONSE  IN  THE 
VICINITY  OF  THE  VERTICAL  AXIS 


To  obtain  a  solution  for  the  radiated  pressure  when  the 
observation  point  falls  on  vertical  axis  <}>  =  0  ,  one  can  substitute 
J^(kKsin,*')  =  1.0  in  Equations  (3.29)  and  (3.31),  yielding: 


I  (R,0, t)  = 


2a5  2a5(3u6iri  -  2a4k2  -  a2k4  +  2a4gr  ,) 


*  exp{P  z  -  ia)  (2  t}  kJ  (k,r)dk 
pi  o  pi  o 


(A.  1) 


1  ( R ,  0 ,  t ) 


ao  S  —  A  ?  A — O  A  _ 

’  /_  _t  ,  2a  Cp2rP2+1[3aX2~2a  k  ~a  k  +2ct  grp2Jl 

.  *•  2o5  2oS(3a 6fl*  -  2a V2  -  aV*  +  2a4gf  >  > 

0  P2 


x  exp { T  z  +  ia)  SI  ,t)  k.T  (kr)dk 
p..  o  pj  o 


(A.  2) 


Using  the  modified  saddle-point  method  and  following  the  same  procedures 

as  outlined  in  Section  3.2.2  to  evaluate  the  above  integrals,  one 

obtains  the  normalized  radiated  pressure  from  an  infinite  elastic  plate 

excited  bv  an  impulse  point  force  Q(t)  =  F  S(t)  given  by: 

o 
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APPENDIX  H 

AP PROX I MATE  SOLUT IONS 


Using  the  methods  of  solution  employed  by  Magrab  and  Reader  [20] 
and  Stuart  [12]  to  evaluate  the  radiated  pressure  from  an  infinite 
elastic  plate  excited  by  a  transient  force,  one  can  compare  such 
solutions  with  the  solutions  obtained  in  Chapter  III.  Applying  the 
inverse  Hankel  transform  to  Equation  (2.20),  one  obtains: 


P(r,r.,w)  = 


r  p  co2Q(U))o*  -'k  J  (kr)dk 

1  O  Q 

V 

/  2  a)  p 

llk  +  -0- 

p 


(B.l) 


2  ii  D 


Letting 


J  (kr)  =  |[H(l) (kr)  +  H(2)(kr) ] 

o  2  o  o 


(B.2) 


H(1)(kr)  x  e 

o  v  krn 


i  (kr- 


(B.3) 


r  =  Rsini]'  ,  z  =  Rcos4>  ,  ft  = 


to 


and  substituting  Equations  (B.2)  and  (B.3)  into  Equation  (B.l)  yields: 


P  (  R ,  4> ,  t ) 


3a4fi2Q(to  -Q)  Jl  e 

_ %  o 

An /hr  si  ni> 


-In 
'  A 


l/9 

k _ “exp{R(rcostj>  +  iksini^)  )dk 

— A  AO  —  '  A  2 

(k  -  a  iT)r  +  Ba  fl 


(B.A) 


—  —2  2  2  1/2 
where  f  =  (k~  -  a  SI  )  ' 


Out*  can  extend  the  real  variable  k  to  the  complex  plane  and  use 


the  saddle-point  method  to  evaluate  the  Integral  In  Equation  (H.4).  The 
choice  ot  the  branch  cuts,  path  of  steepest  descent,  etc.,  are  enumerated 
below. 


Branch  Point  and  Branch  Cuts .  Branch  points  are  defined  by  1'  "  0 
which  yields  k^  =  +  Od!  located  on  the  real  axis.  The  branch  points  and 
branch  cuts  are  shown  in  Figure  B. l  and,  for  detailed  determin.it  ion  ot 
the  branch  cuts,  see  Section  2.1.1a. 

Saddle  Points  and  Path  of  Steepest  Descent  (PSDl  .  The  saddle 
point  (s)  k  of  tin1  integral  in  Equation  (B.'t)  is  defined  as  that  value 
of  k  for  which  the  first  derivative  with  respect  to  k  of  the 
exponential  term  is  zero;  hence, 

k  =  +  otiTsin^  .  (B.5) 

s 

The  saddle  point  and  PSD  are  shown  in  Figure  B.2  and,  for  detailed 
discussion  and  determination  of  PSD  and  poles,  see  Stuart  [12]. 

One  can  use  one  of  two  methods  to  evaluate  these  integrals.  The 
first  method  involves  the  use  of  the  regular  saddle-point  method  which 
does  not  take  into  account  contribution  of  the  pole  when  it  appears  near 
the  saddle  point.  The  second  method  involves  the  use  of  the  modified 
saddle-point  method  which  takes  into  account  contribution  ot  this  pole 
and,  hence,  gives  more  accurate  results  than  the  first. 

(a)  Radiated  Pressure  Obtained  by  Regular  Saddle-Point  Method 
(RSP) .  Evaluating  the  integration  in  Equation  (B.4)  by  the  regular 
saddle-point  method,  one  obtains: 


branch  cuts. 

\  ORIGINAL  PATH 

J INTEGRATION  ^ 

-i<B 

j 

1 

1 

1 

CQ 

_ 

~kl 


Figure  B.  1 
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Applying  the  inverse  Fourier  transform  on  i  irn  in  Fiju.it  ion  (3.6),  Cauchy 
contoux  integration  theorem  is  um  J  to  evaluate  the  integration. 

Figures  B.  )a  atul  H .  5l>  show  poles  and  coniout  of  integration  lor  retarded 
tines  t j  >  0  and  <0  ,  respectively.  It  can  be  shown  that  the 
poles  are: 

E  ,  ■  E.  +  it  ,  E  -  -t  +  it  and  fi  ,»  -i2t~ 

p  1  1  2  p  2  1  J  p  3  J 

where  t  and  s',  are  functions  of  the  angle  s'  • 

The  solutions  for  radiated  pressure  from  an  infinite  elastic 

plate  excited  by  an  impulsive  point  force  Q(t)  =  f  5(t)  become: 

o 

arc  s'  ,F  exp{-2K.u t  /h) 

P(K,^,t)  ■*  - — -  ,  for  t  "*  0  ,  (B.7) 

lTRiriu:  ;si.rv  +  i.o)  1 

ansi  ■■>!/>  K.s;  ,  t  / h 

sloe  (s  +  s~)  F  e  _  „ 

FtK.-.'.t)  •  —  -j  ,  ,  ;  sin{8  +  r  at  s'  -  vl 

2" Rh~  (a'j  +  i\Ty) 

for  T  <  0  ,  ( B . 8 ) 


where 


3(t  j  -  sT,)  sin\‘’ 


2cie: 

arc  t an  (t  , /s  ) 
arc  tan  (a,/a  ) 


T  -  1 


96 


The  normalized  radiated  pressure  is  given  by: 


P(27Th2)  (R/h)  _  2otBoxp [ 2ac^i ^  (R/h)  ] 


Fo(Co/h) 


:  2 - 4 - »  for  T1  1  0  (B . 9) 

[12c2sin  <P  +  1.0] 


P(2xh)  (R/h) 

F  (c  /h) 
o  o 


2  2  1/2  — 

+  e2)  exp[-(R/h)aT]e2] 

7~2~  2.1/2 

(a1  +  a2) 

_ 

x  sin[0  +  -  ax1el  -  i p]  ,  for  x  >  0. 


(B. 10) 


( b )  Radiated  Pressure  Obtained  by  Modified  Saddle-Point  Method 
.(MSPK  Using  Jones'  [21]  modified  saddle-point  method  for  the  evaluation 
of  the  integral  in  Equation  (B.4),  one  obtains: 


_  i3Q(w  ft)^2e 

P(R,t}>,fi)  =  - - 

47T/7TKh 


(aS^sin^  ~  kp)(«ft)  (-2H(Rez)  +  erf c [- (~  ^) 1/2Z] } 
iotficos<j>(ft2sin^<j>  -  1)  +  g 


x  exp  iaSl  +  I(kp  -  cx^sin<]>)2 


where  k  -  pole  close  to  saddle  point  , 

r  1 1/2 


z  =  (k  ^  -  aftsin<(>)  j - — — -■ 

*  ^aQcos2i}i- 


{ — =ii_y  ] . 

OtSlcos“i]r  J 


(B.ll) 


H(x)  =  Heaviside  function. 


Applying  the  inverse  Fourier  transform  to  Equation  (B.ll)  and 


using  Cauchy  contour  integration  theorem,  one  can  evaluate  the  resulting 
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